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Abstract 

In this work we discuss the practical and conceptual issues related 
to quark-hadron-duality in heavy- heavy systems. Recent measure- 
ments in the charmonium region allow a direct test of quark-hadron- 
duality. We present a formula for non-resonant background produc- 
tion in e~^e~ — > DD and extract the resonance parameters of the 
^^(35) — ipi^GS). The obtained results are used to investigate the up- 
silon energy range. 
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1 Introduction 



Quantum Chromo dynamics describes the strong interactions of quarks of 
gluons. Nevertheless, these particles are not experimentally detected as the 
physical states are formed by hadrons. Few methods link the description in 
terms of QCD parameters to the properties of the hadronic bound states. 
Among these are QCD sum rules [HISE], lattice QCD chiral perturba- 
tion theory [OlIZI and the 1/iVc-expansion [Hlini- In this work we focus on the 
method of QCD sum rules where the notion of quark-hadron-duality (QHD) 
plays a dominant role. 

The optical theorem provides the basis for connecting theoretical and 
phenomenological quantities. It relates physical measurable observables like 
the cross section for hadron production to theoretical quantities usually ex- 
pressed by correlators of two- or three-point functions. In the Euclidean 
domain this correlator can be theoretically calculated by means of the op- 
erator product expansion (OPE) JUj. The leading terms are given by the 
perturbative expansion which is supplemented by the condensate contribu- 
tions. This can be compared to the corresponding quantity extracted from 
experiment and in this way it is possible to extract information about the 
system or the QCD parameters. One of the limitations of the sum rules 
already becomes visible. Approaching from a perturbative side, the corre- 
lator does not include real nonperturbative phenomena. Consequently the 
analysis must be performed in a so-called 'sum-rule- window' where the OPE 
of the correlator is under control and the system still reacts sensitive to the 
hadronic parameters. Furthermore, in practical applications the experimen- 
tal spectral density is usually only known for the lowest ground states. To 
estimate the missing information on the phenomenological side, the integral 
over the experimental spectral density is then assumed to equal the inte- 
gral over the theoretical spectral density above a certain threshold energy sq. 
This is the assumption of global quark-hadron-duality. Though being one of 
the basic assumptions in QCD sum rules its range of applicability has only 
been scarcely explored. The foundation of QHD was laid in JT]. Whereas 
in semileptonic decays and lepton scattering the concept of duality is under 
active investigation, e.g. [HlinillllllSlliniini, QHD in the context of QCD 
sum rules has only recently be reinitiated by Shifman, see and references 
therein. 

In this paper we will discuss both the practical and conceptual aspects 
related to QHD. The main part will focus on the charmonium system where 
new measurements from BES ^H] in the region between 3.7 GeV and 4.8 
GeV have improved the experimental situation significantly. Since also the 
theoretical spectral density can be calculated this allows a thorough compar- 
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ison and a stringent test of QHD. We finally extend these investigations to 
the upsilon system. 

The following section is dedicated to an estimate of the uncertainty related 
to the use of QHD. After a discussion of the theoretical contributions we give 
a description for the threshold parameter sq- The error on sq indicates the 
uncertainty related to the assumption of QHD. In section 3 we investigate the 
charmonium cross section in more detail. Apart from the ■0-resonances the 
non-resonant D-production has a significant impact on the cross section. We 
present a model description for this background and extract the resonance 
parameters of the — ip{Q). Section 4 discusses the more conceptual 

issues since the notion of QHD in heavy-heavy-systems is far from trivial. 
The following section concentrates on the upsilon system. In particular, we 
give an estimate for the threshold parameter sq, present a model for the 
non-resonant i?-production and check the validity of the OPE. Finally we 
summarise the results. 

2 Quark-hadron-duality in the charmonium 
system 

In this work we investigate charm production in e"'"e~-collisions 

Rc{s) = ''X'' ^ = 127rlmn(s + ie). (1) 

Via the optical theorem, the experimental cross section is related to the 
imaginary part of the correlator n(s) defined by 

^M) =i\d'x e^'^^ (r{j^(x)jt(0)}) = {q,q. - g,,q') U{q^) . (2) 

The charm vector current is given by in^x) = ec{c'j^c){x) where Cc represents 
the electric charge of the charm quark. 

In principle one can calculate ^{q^) perturbatively, take the imaginary 
part and compare it to the measured cross section. However, a perturba- 
tive calculation of H(g^) is valid only in the Euclidean domain. An analytic 
continuation to the Minkowski region neglects terms which are small in Eu- 
clidean but can become important in Minkowski. Thus, with the assumption 
of global QHD, only smeared quantities can be compared. A further compli- 
cation arises in the theoretical calculation of H(g^). In the deep Euclidean 
domain the perturbative expansion works well. However, usually one is in- 
terested in a region closer to threshold. Here the Coulomb-like behaviour of 
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the charmonium system shows up and the theoretical expansion converges 
badly. These large terms can be resummed with the help of the theory of 
non-relativistic QCD (NRQCD) j2n|l21j. Since we will employ the theoretical 
prediction for n(g^) and the theoretical spectral density in this and the fol- 
lowing chapters, first we briefly discuss these contributions and then return 
to a discussion of QHD. 

The theory of NRQCD provides a consistent framework to treat the prob- 
lem of heavy quark-antiquark production close to threshold. The contribu- 
tions can be described by a nonrelativistic Schrodinger equation and system- 
atically calculated in time-independent perturbation theory. The correlator 
is expressed in terms of a Green's function G{k) = G(0, 0, k) [^1^ 1^: 



where is the number of colours, k = y — s/4 and M^. represents the 
pole mass. The constant Chicts) is a perturbative coefficient needed for the 
matching between the full and the nonrelativistic theory. The contributions 
from NRQCD are summarised in the potential. The Green's function obeys 
the corresponding Schrodinger equation 



Here Vc{x) = —Cpas/l^l represents the Coulomb potential and AV{x) con- 
tains the NLO and NNLO corrections. Details on the solution of this equation 
can be found in [25- The Green's function contains pole contributions be- 
low threshold and a continuum above threshold. In order to construct the 
theoretical continuum spectral density for the full energy range it is not suf- 
ficient to use the spectral density from NRQCD which is only valid for low 
velocities. In addition, one must include the results from perturbation theory 
which gives R{s) at large velocities. At intermediate velocities one can per- 
form a matching between both regimes. In [23] this procedure is described 
in detail. The resulting theoretical spectral density is shown as a solid line 
in fig. ^and we denote the spectral density by R^^'^°{s). 

This spectral density is supposed to give a good approximation to the 
experimental spectral density at high energies. Decreasing the energy, one 
approaches the resonance region. Here R^^^°{s) will fail to reproduce the reso- 
nances. But naively assuming QHD, a smearing of R^^^{s) over a 'sufficiently 
large' energy range should give a good approximation to the same smearing 
with R^^^\s). In fact, the notion of QHD is more subtle and the naive expec- 
tation is not correct. The optical theorem relates the two representations by 




(3) 





(4) 



3 



1 I I I I I I I I I I I I I I I I I I I r 




3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 4.6 4.7 



Figure 1: Cross section from BES with the perturbative spectral density 
(sohd) and the background contribution (dashed). 

a dispersion relation which includes all values of the energy. So only the com- 
plete phenomenological result can be compared to the full theoretical result 
which also includes the poles of the Green's function. In section 4 we discuss 
the concept of QHD in more detail and present quantitative estimates for 
the individual contributions. In fig. ^we have plotted and the mea- 

sured cross section in the energy range between 3.7 GeV < ^/s < 4.8 GeV. 
At these energies clearly lies above the data points. On the other 

hand, the contribution from the theoretical poles turns out to be smaller 
than the contribution from the lowest ■i/'-states. 

To test the accuracy of QHD from a comparison of the theoretical and 
phenomenological spectral densities we choose the moments 

Wn{A,so,s,) = r ds , ^^^} . (5) 

The weight function w{s) = l/{s + A)""*"^ corresponds to the one usually 
used in the moment sum rules to extract the quark masses. Another popular 
weight function is Wu{s) = e~^/" which is used in the Borel sum rules. Since 
the analysis and the results are very similar in both cases we will not perform 
an independent analysis for the Borel sum rules in this work. A and n 
are free parameters which can be used to move the sum rules to a region 
convenient for the analysis. Large A will improve the perturbative expansion 



4 



while for small or even negative A the analysis will react very sensitive to 
the bound states. Small n will result in a relative fiat weight function and 
high n put the emphasis on the low energy region. As will be discussed in 
chapter 4, si should be taken to infinity to establish a precise relation between 
the phenomenological and theoretical part where they are connected by an 
Euclidean quantity. However, in practice the cross section is only measured 
up to a certain energy so above this energy one has to rely on the theoretical 
prediction for R{s). When comparing both parts the integral above this 
energy is then equal. 

The experimental cross section can be interpreted as non-resonant back;- 
ground production of D'^^^-mesons and resonances of Breit-Wigner form. For 
a comparison to R*^'^° it is interesting to separate these two contributions. 
In the next section we will give a model description of the background pro- 
duction which has been plotted as a dashed line in fig. [T] and discuss the 
charmonium cross section in more detail. In this section we focus our at- 
tention on the most important issue in actual sum rule calculations and the 
basic question of QHD: how well can the experimental cross section be ap- 
proximated by the theoretical spectral density and how large is the error on 
iyn(A, So, Si)? The BES data have been given directly in terms of B^^^{s) 
and we will use these values as our reference for the experimental cross sec- 
tion. From the measured spectral density the light quark contribution must 
be subtracted. At these energies the light quarks can safely assumed to be 
massless and the high energy approximation [211 provides a good descrip- 
tion. Apart from BES, the resonance properties of the V'-states have been 
extracted in [2Z1I2H! and of the T-states in [23120! • 

In order to test a typical problem in the sum rules let us now assume that 
the only information we had from the experimental side were the resonance 
properties of the first p bound states. With the assumption of QHD, the 
contribution from the higher states is then given by the integration of the 
theoretical spectral density above a threshold sq. The above question 

can then be formulated in a different way: what value of sq must be taken? 

Since the QCD sum rules are used to extract the heavy quark masses, 
the choice of Sq influences the central value of the masses. Furthermore 
the uncertainty in sq translates directly to the error of the masses and it is 
therefore important to have a reliable estimate of this uncertainty. There 
is no rigorous justification for a particular choice of sq. As a heuristic rule 
it is usually assumed that Sq should be given by about 250 MeV above the 
highest included resonance. In this section we want to check if - and to what 
extend - this rule is valid. It will turn out that particular care must given 
to a possible background contribution. In chapter 5 we will then apply the 
results obtained in this section to estimate sq in the upsilon system where 
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3.097 


3.686 


3.770 


4.040 


4.159 


4.415 




n 


= 


3.78 


4.07 


4.10 


4.21 


4.31 


4.38 


v^/[GeV] 


n 


= 3 


3.65 


3.99 


4.03 


4.14 


4.25 


4.30 
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= 7 


3.59 


3.92 


3.99 


4.08 


4.17 


4.20 
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= (4 GeVy 
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5 


6 




n 


= 


3.82 


4.09 


4.12 


4.22 


4.32 


4.39 


v^/[GeV] 


n 


= 3 


3.72 


4.03 


4.07 


4.18 


4.29 


4.35 




n 


= 7 


3.65 


3.98 


4.03 


4.13 


4.24 


4.30 



Table 1: sq for different number of poles p and different values of n and A 
from eq. (jHI). 



no experimental information is available to fix this parameter. 

To determine sq in the charmonium system we now compare the quantities 

fc=i (El + A)-+i ^ 73.72 (s + A)"+i 

h (El + A)"+i Js[r^ {s + A)"+i ' 

where a is the electromagnetic fine structure constant, Tee,k is the partial 
decay width into e^e~ and Ek is the mass of the k*^^ resonance. We have 
used the narrow-width approximation for the resonances. The sum over the 
resonances in the first line extends only over the first 2 resonances since the 
other known resonances ip{3) — ip{6) are included in R^^^{s). The upper 
limit of the integration is taken to infinity, si oo. Since the BES data 
have only been measured up to 4.8 GeV we assume QHD above this energy 
so the integral from 4.8 GeV to infinity is identical on both sides and drops 
out. The first line represents the 'exact' result from the data. The second 
line is a typical phenomenological approximation using the assumptions of 
QHD in which sq is left as a free parameter. As mentioned above, it is usually 
assumed that sq should be given by the mass of the highest resonance Ep 
plus about 250 MeV. In table [T] we have listed Sq obtained from eq. (0) for 
different values of n and A. 

Let us first look on the behaviour of Sq on n and A. We see that Sq depends 
on the choice of these parameters. Since R^^'^" lies above the experimental 
cross section, larger values of n will lower sq. The difference between the 
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largest {n = 0) and smallest {n = 7) value varies between A^/sq = 190 MeV 
for p = 1 and = 110 MeV for p = 3 at A = 0. At A = (4 GeV)^ the 

analysis is in a more perturbative region. Therefore one expects less impact 
of n on So, but still the difference remains sizeable: from A^/sq = 170 MeV 
for p = 1 to A^/so = 80 MeV for p = 5. This is a remarkable result for 
applications of the sum rules. In analyses where the threshold Sq has an 
important impact on the quantity one would like to extract, the change of 
So with n and A might influence the final result. 

Furthermore we note that the rule 'highest resonance plus 250 MeV is 
strongly violated. Taking only the lowest pole, p = 1, this is no surprise. 
Since the first two poles are very dominant on the experimental side, one 
cannot hope to give a good description of these poles by the perturbative 
spectral density without taking into account the pole contributions from the 
Green's function. Also the reason for the violation of the rule for the higher 
states is clear: Using only the resonance parameters Ek and Tee,k in the 
second line of eq. © we have neglected the non-resonant D*^*^-production. 
There are two ways to estimate this background contribution. The first one 
is to give a model description for the background as has been depicted in fig. 
m In a second approach one could assume that QHD above the third pole 
already represents a reasonable description. In this case the phenomenolog- 
ical part is given by '3 poles plus R^^'^° from a threshold of 250 MeV above 
the 3'"'^ pole'. Without background this description also applies to p poles 
instead of three. Subtracting these two descriptions should therefore give an 
estimate of the background contribution. 

The drawback of both methods is clear: in the first one the result is model 
dependent where in the second one it was assumed that QHD could already 
be used for states with p > 3. To include the background contribution 
one could either add this background explicitly to the second part of eq. 
(jHI) and use ^/so = Ep + 250 MeV or stick to eq. (jHl) and lower sq by 
the appropriate value. Since we want to compare the results with table 
H we use the second method, sq is then determined from the equation 

dsR^^^°{s)/{s + A)"+^ = 'background' where ^ = Ep + 250 MeV and 
the estimate of the background ranges up to Si. The results for so are shown 
in table |21 In the first row the background has been estimated by the model 
description from fig. Q and the second row in parentheses shows the result 
assuming QHD already for p > 3. 

Let us now compare the 'exact' result for so from tabled which is needed 
for a correct description of the experimental moments to the QHD-based 
estimate of So in table El For p = 1 QHD cannot reproduce the correct value. 
This is expected since the theoretical spectral density will not give a good 
description for the second resonance. At p = 2, 3 the background has almost 



7 



A = 
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n = 





3.35 


3.93 


3.99 


4.16 


4.22 


4.32 








(3.35) 


(3.94) 


(4.02) 


(4.12) 


(4.22) 


(4.29) 


x/5^/[GeV] 


n = 


3 


3.35 


3.93 


3.99 


4.14 


4.18 


4.25 








(3.35) 


(3.94) 


(4.02) 


(4.12) 


(4.23) 


(4.28) 




n = 


7 


3.35 


3.93 


3.99 


4.11 


4.14 


4.17 








(3.35) 


(3.94) 


(4.02) 


(4.13) 


(4.24) 


(4.28) 


A = (4 GeV)2 
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2 
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4 


5 


6 




n = 





3.35 


3.93 


3.99 


4.16 


4.22 


4.33 








(3.35) 


(3.94) 


(4.02) 


(4.12) 


(4.22) 


(4.29) 


v^/[GeV] 


n = 


3 


3.35 


3.93 


3.99 


4.15 


4.21 


4.29 








(3.35) 


(3.94) 


(4.02) 


(4.12) 


(4.22) 


(4.28) 




n = 


7 


3.35 


3.93 


3.99 


4.14 


4.18 


4.24 








(3.35) 


(3.94) 


(4.02) 


(4.12) 


(4.23) 


(4.28) 



Table 2: sq for different number of poles p and different values of n and A 
obtained from the description 'a/sq = Ep plus 250 MeV minus background'. 
The numbers in the first row are with a background from fig. C]and the num- 
bers in parentheses with the second estimate of the background as described 
in the text. 
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no influence on sq so it essentially lies 250 MeV above the resonance. For 
a relatively steep weight function, n = 7, the results are similar to table [H 
but for n = they differ up to 160 MeV. The estimates for sq agree well 
for p = 4 and for p = 5,6 they are a bit lower. In general, we see that the 
assumption of QHD with the corresponding choice for sq gives a reasonable 
description of the experimental moments, at least for the higher poles. It is 
interesting to note that in table |21 the change of Sq with n and A is relatively 
small. With the background model it is smaller than in tabled especially for 
a small number of poles p. In the second approach with QHD assumption for 
p > 3 the value of sq remains almost constant. So this variation can easily 
be underestimated. 

To summarise, we have estimated the uncertainty connected with the 
use of QHD. In the charmonium system QHD represents a reasonable good 
approximation if at least the first two poles are added explicitly to the phe- 
nomenological side. The moments can then be determined with the de- 
scription 'poles with resonance parameters plus background plus theoretical 
spectral density above Sq'. The threshold Sq should be given by the energy of 
the highest pole Ep plus 250-300 MeV. Estimating the uncertainty on sq one 
should take into account that the variation with n and A can easily amount 
to 100 MeV. In addition one should allow a variation of 100 MeV from its 
value ^Ep plus 250-300 MeV. So we conclude that a reasonable error estimate 
for So is given by A^/sq = 200 MeV around its central value. 

3 Charmonium cross section 

In this section we investigate the charmonium cross section in more detail. 
The first two resonances are dominating the cross section clearly, the J/ip 
at 3.097 GeV with a partial decay width of Fee = 5.26 keV and the ip' at 
3.686 GeV with Fee = 2.19 keV. Since both resonances lie below open D- 
production their total widths are small, Ttot = 87 keV and Ttot = 300 keV [HI] 
(264 keV [32J) respectively. At 3.74 GeV open Z^Z^-production starts. In 
the continuum 4 more resonances have been identified: a relatively small 
resonance ip{3.770) just above Z)D-threshold and three broader resonances, 
■?/'(4.040), ijj{4.159) and '?/'(4.415). From the data it can be seen that in the 
energy range above 4 GeV the background continuum gives a significant 
contribution. 

Now we want to give a model description for the background. Our mo- 
tivation is twofold. As already seen in the last section, it is interesting to 
separate the background and resonance contributions to estimate the relative 
size and the importance of the higher resonances. Furthermore, in practi- 
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cal applications it is more convenient to deal with a smooth approximating 
function in terms of a few resonance parameters than with a large number 
of data points. 

The next channel above DD starts at 3.88 GeV with DD and D*D- 
pro duct ion. At higher energies open the DsDs, D*D , D^Ds , and D^Ds 
channels. Since D* has three spin directions, the production of DD and 
D*D is enhanced compared to DD-production • 

One could try to parametrise all these contributions by appropriate form 
factors. However, neither theory nor experiment provide sufficient informa- 
tion to predict these form factors. As a consequence, this ansatz would 
depend on many free parameters that had to be fitted from the used data. 
So the result would strongly depend on the data set and could not be gen- 
eralised. Therefore we use a different approach based on perturbative QCD 
and model the non-resonant background production by 



For as{s) we use the three-loop formula with A = 313 ± 27 MeV which 
corresponds to as{Mz) = 0.1172 ± 0.002 and three light flavours. This 
background is plotted as a dashed line in flg. ^ For energies sufficiently above 
cc-threshold we expect that the main process will be hard cc-production 
which finally turns into D-mesons with unit probability. The higher order 
corrections in eq. (jZj) are based on perturbation theory in the massless limit 
j2ni and show the correct high energy behaviour. For a finite charm quark 
mass the expansion contains soft gluons which are interchanged between the 
quarks at threshold. These soft gluon ladders, which can be resummed by 
means of NRQCD, lead to the formation of the ■j/'-states and are responsible 
for the resonance effects. Thus they must not be included in the background 
description. 

It remains the choice of the threshold Mfhre which appears in the phase 
space factor. In a QCD-based picture it is given by the charm pole mass 
Mthre = 2Mc. However, in this case we describe Z)-meson production, so the 
phase space should rather be given by the phase space the D-meson than that 
of the charm mass. As described above, the production of DD is suppressed 
to DD . In addition, for higher energies the other heavier thresholds open. 
Consequently a threshold of Mthre = M^^^ would probably overestimate the 




Mthre 



rDD* 

thre 



3.879 GeV. 



(7) 
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k 


Ek/[GeV] 


rL/[MeV] 


rL/[keV] 


3 


3.7727 ±0.0016 


24.4 ±4.3 


0.19 ±0.025 


4 


4.0504 ±0.0043 


98.5 ± 12.8 


1.03 ±0.11 


5 


4.1665 ±0.0061 


55.9 ± 12.3 


0.37 ±0.081 


6 


4.4294 ±0.0085 


86.0 ±20.9 


0.39 ± 0.074 



Table 3: E^, F^^t and Fgg for the resonances ipiS) — '?/'(6) from eq. (jHI) and 
using BES data ^ between 3.7 GeV and 4.6 GeV. 



background in the intermediate energy range. Therefore we fix the threshold 
parameter to Mthre = M^^J = 3.879 GeV. One of the limitations of this 
description is obvious: one cannot expect it to be a good approximation 
below and directly above the DD -threshold in the range of 3.74 GeV < ^/s < 
4.0 GeV. Below 3.88 GeV we obviously miss the small DZ)-production which 
may also have some impact directly above the DD threshold. However, we 
expect that ((Zj) gives a good description of the background in the energy 
region between 4.0 GeV < a/s < 4.6 GeV. Above y/s > 4.6 GeV the 
resonance structure seems to level off into a continuum. Here one cannot 
separate the background and resonances any longer. Thus, instead of the 
background, one should use the full theoretical result to describe the spectral 
density for energies above 4.6 GeV. 

Now we give a description of the cross section in terms of the background 
and Breit-Wigner resonances: 

We prefer to use a constant total width F^^j instead of a s-dependent one 
since the functional form close to threshold is not clear and the uncertainty 
connected with the background description is at least of the same order. 
With this formula one can extract the resonance parameters from the BES 
data ^H! between 3.7 GeV and 4.6 GeV (75 data points). The statistical and 
systematic error have been added quadratically for each data point. In table 
Elwe have hsted our results for E^, F^^^ and F^^ for the resonances '?/'(3)— ■?/'(6). 
The fit gives a x^/ dof = 1.02. In fig. |2lwe have plotted the resulting spectral 
density together with the data points. It can be seen that the experimental 
spectral density is rather well approximated by the theoretical description. 
We can compare the results to former measurements where the estimate of 
the non-resonant D-production was fitted to the data j2Zll2HlEIl- We see 
that the masses change only by several MeV. Whereas the change of F^^^ and 
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Figure 2: Spectral density from (jHl) with resonance parameters from tab. IHl 
and BBS data HH]. 

for the '?/'(3) is mild, these changes are larger for the higher resonances 
since these parameters are not only based on a different data set but also on 
a different background description. 

To investigate closer the dependence of the resonance parameters on the 
background we now vary the threshold energy Mt^re- We compare to the 
results obtained from a background with Mfhre = 3.74 GeV which would 
correspond to Dl)-production and to Mthre = 4.0 GeV. Compared to the 
data, the first one seems to overestimate and the second one to underestimate 
the background. 

We have listed the results in table 01 The statistical errors are very 
similar to table El and have therefore been omitted. The masses remain very 
stable. The masses of the ipiS) and ipi^) do not change, ipi^) changes by 
Ai?5 = ±3.5 MeV and the largest change is of the ■?/;(4) with AE4 = -5 MeV 
for Mthre = 4.0 GeV. For Mthre = 3.74 GeV the widths of the tp{3) are 
reduced since the background already starts below the resonance. The widths 
of the seem to be very stable against variation of the background and 
the change for the ip{Q) is of a similar size as the statistical errors. The 
widths of the ip{4:) are the most sensitive to a variation of the background. 
They show a significant change which is clearly larger than the statistical 
error. These changes of the resonance parameters are part of the systematic 
uncertainty connected with the background description and of the difficulty 
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Ek/[GeV] 


tot/ 


MeV] 


eel 


keV] 


k 


a 


b 


a 


b 


a 


b 


3 


3.7715 


3.7725 


18.0 


22.4 


0.12 


0.17 


4 


4.0504 


4.0450 


75.5 


138.1 


0.74 


1.60 


5 


4.1629 


4.1701 


56.4 


51.5 


0.36 


0.32 


6 


4.4296 


4.4299 


70.8 


98.3 


0.30 


0.48 



Table 4: Resonance parameters for different backgrounds with M^/ire = 
3.74 GeV (a) and Mthre = 4.0 GeV (b). The statistical errors are similar to 
table El 



to separate these two contributions. However, the experimental cross section 
is well approximated by eqs. ()7p8|l and the resonance parameters of table 01 

4 Theoretical versus phenomenological mo- 
ments 

QCD sum rules provide a framework which relates a QCD-based description 
in terms of QCD parameters to measurable quantities in terms of hadron 
properties. In this section we discuss the different conceptions related to the 
theoretical and phenomenological description. 

The advantage of using a weight function as in ^ is that the moments 
Wn are directly connected with n(s) at the Euclidean point s = —A: 

127r^ 

WniA)=Wn{A,Srmn,00) = ^—n{s) (9) 

n! as" s=-A 

where Smin indicates the lowest pole. In the Euclidean region the theoretical 
expansion is known to be valid by means of the OPE. In addition to the 
perturbative result, condensates of higher and higher power will appear. In 
a QCD-based picture the definition of the pole mass provides the natural 
description for the onset of the continuous spectral density and therefore 
the threshold is given by y/s = 2M(.. The theoretical expansion depends on 
the values of A and n: large values of A and small n move the moments 
to a safe perturbative region and the expansion in as converges well. In 
principle, Wn could thus be calculated to high accuracy. However, this region 
is of little phenomenological interest. Usually in sum rules analyses one is 
interested in extracting information on the ground state or the quark masses. 
In order to be sensitive to these parameters the analysis must be performed 
relatively close to threshold. In this case the perturbative expansion does 
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not converge well any more since large terms appear reflecting the Coulombic 
structure of the charmonium system. These potentially large terms can be 
resummed with the method of NRQCD which sets up a systematic framework 
to treat these non-relativistic corrections. The result is expressed in terms of 
a Green's function and can be directly evaluated at s = —A. Its imaginary 
part shows poles below and a continuum spectral density above threshold. 

This QCD-based theoretical description has to be confronted to the mea- 
sured cross section to which it is related by the optical theorem. As described 
in the last section, its behaviour is very complicated: It contains two sharp 
resonances and several higher resonances which are shifted into the contin- 
uum. It is obvious that the theoretical and phenomenological spectral density 
do not equal each other. In the Euclidean region the theoretical expansion is 
truncated - as a series in as and in higher condensates. If one were able to 
calculate n(g^) in the Euclidean domain exactly, one could analytically con- 
tinue the result to the Minkowski domain and take the imaginary part. The 
theoretical spectral density would equal the hadronic cross section. However, 
in practice only the truncated expansion is analytically continued to values 
of positive q^. Small neglected terms in the Euclidean domain can become 
large in the Minkowski region and change the spectral density significantly. 
The origin and behaviour of such contributions have been discussed in |18j . 

In the OPE condensates appear, the leading contribution is given by the 
gluon condensate: 



the analytic form of the functions Cpp and Cpp can be found in |^. How- 
ever, its contribution to the moments is small for values of n and A used 
in this analysis. Its contribution grows if one comes very close to thresh- 
old, for A > — 4M^ or for very large n since here the moments test the 
nonperturbative region. 

Now we want to compare the size of the individual theoretical and phe- 
nomenological contributions. In fact, this comparison is done in QCD sum 
rules to extract the charm and bottom quark masses since the moments 
show a strong dependence on the value of the mass. We fix the MS-mass 
to mc(mc) = 1.19 GeV j25j. For the comparison we use a range of val- 
ues for n and A somewhat larger than in typical sum rule applications. In 
table El we show the results for the theoretical and experimental moments. 
The moments have been normalised to the total phenomenological moments. 
The poles represent the dominant contribution and they are even more pro- 
nounced on the phenomenological side. It can be clearly seen that small A 



4 {^FF) 



{Cfp{s) + y?pis)) , 



(10) 



9 16M4 
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A 


= 






Ti 


3 


7 


12 


Theory total 


1.04 


0.85 


0.67 


Theory poles 


0.71 


0.76 


0.65 


Theory continuum 


0.33 


0.086 


0.021 


Exp. total 


1 


1 


1 


Exp. poles 1+2 


0.87 


0.989 


0.9992 


Exp. BES data 


0.13 


0.011 


0.0008 


A = (4 GeV)2 


n 


3 


7 


12 


Theory total 


0.98 


1.10 


0.99 


Theory poles 


0.34 


0.72 


0.80 


Theory continuum 


0.64 


0.38 


0.19 


Exp. total 


1 


1 


1 


Exp. poles 1+2 


0.56 


0.86 


0.96 


Exp. BES data 


0.44 


0.14 


0.04 



Table 5: Charmonium: Theoretical and phenomenological moments nor- 
malised to the total phenomenological moments. 



and large n shift the analysis closer to the poles. For n = 12 the continuum 
region is essentially cut off. For this value of the charm mass the moments 
show a good stability, which is no surprise since this stability criterion was 
used to extract the value of the charm mass The total theoretical mo- 
ments differ from the phenomenological ones by about 10%. The convergence 
is better for large values of A. Only for A = and n = 12 they differ sig- 
nificantly, but here the theoretical moments are evaluated close to threshold 
and one cannot expect a reliable description of these moments. 

It is interesting to investigate the dependence of the theoretical moments 
on the mass. With a MS-mass of mc(mc) = 1.3 GeV the total theoretical 
moments normalised to the phenomenological ones for A = vary from 0.82 
(for n = 3) to 0.21 {n = 12) and from 0.96 {n = 3) to 0.68 {n = 12) at 
A = (4 GeV)^. For mc{mc) = 1.1 GeV the relative theoretical moments for 
A = increase from 1.33 {n = 3) to 2.27 [n = 12) and from 0.97 [n = 3) 
to 1.4 {n = 12) at A = (4 GeV)^. It can be seen that the sensitivity of the 
moments on the mass decreases for large values of A. The change is due 
to the pole contribution since the continuum spectral density is practically 
independent of the mass. These results confirm the significant effect of the 
mass on the moments and on the stability of the sum rules. 
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5 Upsilon system 



In the upsilon system the experimental situation is unsatisfactory. Apart 
from the resonance parameters of the 6 T-states [31J almost no direct infor- 
mation on the cross section above the T(4) is available. Therefore it is not 
possible to perform a check of QHD as has been done in section 2. However, 
we can use the results of the last sections to draw some conclusions for this 
energy region. 

First we compare the size of resonance and background contribution. 
Analog to eq. (|7j) we make the following ansatz: 

X f 1 + ^ + 1.4092 ( ^] ' - 12.805 ( ^ 

Mthre = M,^f; = 10.604 GeV. (11) 

The strong coupling constant at three loop with 4 light flavours is determined 
from A = 278 ± 30 MeV. Similar to the charm case, we fix the start of the 
threshold at the energy for i?i?*-production. 

We can compare the background to the resonance contribution. In 
fig. El we have plotted the background versus the spectral density using the 
resonance parameters T(4) — T(6) from ^T]. The data indicate P^llHUj that in 
the region between 11.0 GeV < ^/so < 11.2 GeV the cross section continues 
at values around the peak of the T(6). It can be seen that the resonance 
contribution exceeds the background contribution in the range until ~ 11.05 
GeV. In the resonance region the influence of the background is not as strong 
as in the charmonium system. Above 11.05 GeV, since no further resonances 
are measured, one should use the full theoretical spectral density to describe 
the cross section. 

An important parameter in sum rule calculations of the bottom mass is 
the threshold parameter sq. The phenomenological moments are determined 
from the first six T-resonances and the theoretical spectral density above 
Sq. With no background production this parameter would be estimated to 
lie about 250 MeV above the T(6). However, the non-resonant background, 
where it is not already included in the resonance parameters, will effectively 
lower this threshold. In chapter 2 two ways to estimate this parameter have 
been discussed, by a model description for the background and by assuming 
QHD already for states n > 3. In the first one we use the background of 
eq. (fTT|). In the region between 10.75 — 11.05 GeV the total cross section is 
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Figure 3: T(4) — T(6) with theoretical spectral density (solid) and back- 
ground (dashed). 



A = 


= 






n 


3 


7 


12 


Theory total 


1.05 


0.94 


0.76 


Theory poles 


0.41 


0.60 


0.61 


Theory eoiitiiiuuiii 


0.61 


0.31 


0.15 


Exp. total 


1 


1 


1 


Exp. poles 1+2 


0.70 


0.94 


0.991 


Exp. BES data 


0.30 


0.06 


0.009 


A = (10 GeVf 


n 


3 


7 


12 


Theory total 


1.04 


1.06 


0.99 


Theory poles 


0.22 


0.45 


0.58 


Theory continuum 


0.82 


0.61 


0.41 


Exp. total 


1 


1 


1 


Exp. poles 1+2 


0.45 


0.75 


0.91 


Exp. BES data 


0.55 


0.25 


0.09 



Table 6: Upsilon: Theoretical and phenomenological moments normalised 
to the total phenomenological moments. 
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well approximated by the 5 and 6 resonance. Calculating the background 
for the remaining interval the effective threshold lowers to ^/sq = 11.15 — 
11.20 GeV. In the second approach we assume QHD for states with n > 3. 
The background from this ansatz would reduce the threshold to ^/so = 11.0 — 
11.05 GeV. Therefore, taking the average between both estimates, we suggest 
that a reasonable estimate for sq would be given by ^/so = 11.1 ± 0.2 GeV. 

As in section 4 we can compare the contributions from the theoretical 
and phenomenological moments. Again we evaluate the moments in a re- 
gion where threshold effects are important. In table IHl we have collected the 
moments from the different sources. As input parameter for the mass we 
have used mh{mb) = 4.24 GeV |2S]. For the employed values of n and A the 
dominance of the poles is less pronounced than in the charmonium. The the- 
oretical moments differ by about 5% from the phenomenological ones. Thus 
they show a better convergence than in the charmonium system. At A = 
and n = 12 the theoretical calculation of the moments should not be trusted 
any more. 

Now we look at the influence of the bottom mass on the moments. With 
a MS-mass of mi,{mh) = 4.34 GeV the total theoretical moments for A = 
vary from 0.93 (for n = 3) to 0.47 {n = 12) and from 0.99 {n = 3) to 0.80 
{n = 12) at A = (10 GeV)^. For rribimb) = 4.14 GeV the moments for A = 
increase from 1.19 {n = 3) to 1.23 (n = 12) and from 1.1 {n = 3) to 1.24 
(n = 12) at A = (10 GeV)^ thus showing again the strong influence of the 
quark mass. 

6 Conclusions 

In this work we have explored the cross section of e"'"e~-collisions in the 
charmonium and upsilon energy region. The experimental situation for the 
charmonium has significantly improved with the new results from BES [TH] . 
This allows a thorough comparison of the theoretical and experimental de- 
scription for this energy range. These investigations have been applied to 
heavy-heavy systems. Apart from the general discussions it is not clear if the 
concrete results could be generalised to e.g. heavy-light or light-light states 
since the underlying physical systems are different and their properties are 
rather determined from the dynamics of the light quarks. 

The main part of the paper has investigated the charmonium energy 
range. We have given a prescription for the threshold parameter sq which is 
needed to describe the experimental spectral density for large energies. The 
error on this quantity was estimated to Ay/so = 200 MeV which includes a 
possible variation of sq with n and A. The phenomenological cross section 
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can be described as Breit-Wigner resonances and a non-resonant background 
production of D(*)-mesons. We have presented a model for this background 
production based on perturbative QCD. The masses, the hadronic widths and 
the partial e"'"e~-widths of the states iIj{3)—iIj{6) have been extracted from the 
BES data. So we obtain a direct theoretical description of the experimental 
cross section in terms of few resonance parameters with a x^/dof = 1.02. 
QCD sum rules use the identity of the theoretical and phenomenological 
moments related by the optical theorem and can be used to extract e.g. 
the quark masses or the ground state properties. In order to be sensitive to 
these parameters the moments must be evaluated in a region not too far from 
threshold. Thus the contributions from NRQCD form an essential part on the 
theoretical side reflecting the fact that the underlying system is a Coulombic 
one. The different contributions from the poles and the continuum part have 
been compared for different values of n and A. 

Section 5 has been devoted to the upsilon system. Here the experimental 
situation is dissatisfactory. However, the results of the previous sections could 
be used to investigate several properties of this system: a model description 
of non-resonant S(*)-production has been presented and the threshold pa- 
rameter So estimated. We have compared the different contributions to the 
theoretical and phenomenological moments and investigated the effect of the 
bottom mass. Unfortunately the cross section is not well measured in the 
region above the T(4). A more detailed knowledge would allow a better test 
of QHD in this energy region. 

The two basic pictures of QCD are related by QHD: the hadronic world 
and a description based on perturbative QCD in terms of quarks and gluons. 
Thus a better understanding of QHD could provide further insight into the 
structure and behaviour of nonperturbative contributions. 
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